AN EXPLICIT FINITE DIFFERENCE SCHEME 
FOR THE CAMASSA-HOLM EQUATION 

G. M. COCLITE, K. H. KARLSEN, AND N. H. RISEBRO 

Abstract. We put forward and analyze an explicit finite difference scheme for 
the Camassa-Holm shallow water equation that can handle general initial 
data and thus peakon-antipeakon interactions. Assuming a specified condition 
restricting the time step in terms of the spatial discretization parameter, we 
prove that the difference scheme converges strongly in towards a dissipative 
weak solution of Camassa-Holm equation. 
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1. Introduction 

In this paper we present and analyze an explicit finite difference scheme for the 
Camassa-Holm partial differential equation [7] 

(f.l) dtu - df^^u + 3uda:u = 2d,,udl^u + ud^^^u, {t,x) e {0,T) xR, 

which we augment with an initial condition: 

(f .2) u\t=Q = ua, ua £ H^R), ua ^ 0. 

Rewriting equation (|l.ip as 

(I - dlJ [dtu + ud^u] + a, (u^ + ^{d,u)A = 0, 
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we see that (for smooth solutions) (|l.ip is equivalent to the elliptic-hyperbolic 
system 

(1.3) d,u + ud.u + d.P = 0, -aLP + P = + lid^uf. 

Recalling that e~l^l/2 is the Green's function of the operator 1 — d^^, (|1.3|) can be 
written as 

(1.4) dtu + d^F{u,d,u) = 0, F{u,d,u) = ^ 

which can be viewed as a conservation law with nonlocal flux function. In this paper 
the relevant formulation of the Camassa-Holm equation (jl.ip is the one provided 
by the hyperbolic-elhptic system (|1.3|) or (|1.4[) . 

The Camassa-Holm equation can be viewed as a model for the propagation of 
unidirectional shallow water waves [71 [32] ; it is a member of the class of weakly 
nonlinear and weakly dispersive shallow water models, a class which already con- 
tains the Korteweg-de Vries (KdV) and Benjamin-Bona-Mahony (BBM) equations. 
In another interpretation the Camassa-Holm equation models finite length, small- 
amplitude radial deformation waves in cylindrical compressible hyperelastic rods 
[21] . It arises also in the context of differential geometry as an equation for geodesies 
of the if ^-metric on the diffeomorphism group, see for example [171 [HI [30l [36] . The 
Camassa-Holm equation possesses several striking properties such as an inifinite 
number of conserved integrals, a bi-Hamiltonian structure, and complete integra- 
bility in [71 [191 [SI [IS]- Moreover, it enjoys an infinite number of non-smooth 
solitary wave solutions, called peakons, which are weak solutions of p.4|l . 

From a mathematical point of view the Camassa-Holm equation has by now 
become rather well-studied. While it is impossible to give a complete overview of 
the mathematical literature, we shall here mention a few typical results, starting 
with the local(-in-time) existence results in [T5l[34l[37] and those using Besov spaces 
in [23l [22] . It is well-known that global solutions do not exist and wave-breaking 
occurs [7]. Wave-breaking means that the solution itself stays bounded while the 
spatial derivative becomes unbounded in finite time. 

In view of what we have said so far (peakon solutions/ wave-breaking) it is clear 
that a theory based on weak solutions is essential. In the literature there are a 
number of results on (dissipative and conservative) weak solutions of the Camassa- 
Holm equation, see [31 [H [3 [Ml El |2Ql [29l [39l [40] and the references cited therein. 
In this paper we are interested specifically in the class of dissipative weak solutions 
studied by Xin and Zhang [321 HO]- Their results show, among other things, that 
there exists a global dissipative weak solution of (|l.ip - (|1.2p for any initial data 
uo (peakon- antipeakon interactions are covered). These solutions are global in the 
sense that they are defined past the blow-up time (wave- breaking) . More precisely, 
suppose Mo G H^{M.). Then there exists a global weak (distributional) solution 
u e L°°(0,T;7Ji(M)) of satisfying the following properties: t ^ \\u(t, •)|liji(K) 
is non-increasing; dxU G ifQ^(M+ x M), p < 3; 

2 

(1.5) i9:ru(t,a;) < - + C||uo||jfi(R) , for <> 0, 

for some positive constant C. This last item presumably singles out a unique weak 
solution. As an example of how this may work we consider the "peakon-antipeakon" 
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solution given by 

(1.6) u{t, x) = tanh(i - 1) (^5-^^-^^'-^^^ - e-l^+2'(*-i)l) , y{t) = log (cosh(i)) . 

This formula represents a peakon (e"'^"''^!) colliding with an antipeakon (— e^'^^^l) 
at a; = and t = 1. Note that u{l, x) — 0. How this solution is extended to t > 1 
depends on which solution concept we adopt. If we use the formula (|1.6p also for 
t > 1 we get the conservative solution for which \\u(t, •)II//i(r) is constant for almost 
all t. We can also extend the solution by defining u{t, x) = ioi t > 1. Obviously, 
the "entropy condition" ()1.5p will only be satisfied for this dissipative solution. 

Let us now turn to the topic of the present paper, which is the design and analysis 
of numerical schemes. The first numerical results for the Camassa-Holm equation 
are presented in [8] using a pseudo-spectral scheme. Numerical simulations with 
pesudo-spectral schemes are also reported in [25l [31] . Numerical schemes based on 
multipeakons (thereby exploiting the Hamiltonian structure of the Camassa-Holm 
equation) are examined in [SI [SI [TO]- In [lH]j the authors prove that the multipeakon 
algorithm from [21 [TO] converges to the solution of the Camassa-Holm equation as 
the number of peakons tends to infinity. This convergence result applies to the 
specific situation where the initial function uq G is such that (1 — d^,j.)uQ is 
a positive measure. For the same class of initial data, in the authors prove 
that a semi-discrete finite difference scheme based on the variable m := (1 — d^^)u 
converges strongly in to the weak solution identified in [TO] [2D]. In [S^, the 
authors establish error estimates for a spectral projection scheme for smooth solu- 
tions. In a different direction, an adaptive high-resolution finite volume scheme is 
developed and used in [1]. The local discontinuous Galerkin method is adapted to 
the Camassa-Holm equation in [41] . Although this work does not provide a rigorous 
convergence result for general (non-smooth) solutions, they show that the discrete 
total energy is nonincreasing in time, thereby suggesting that the approximate 
solutions arc of dissipative nature. Besides, they establish an error estimate for 
smooth solutions. Finally, multi-symplectic schemes possessing good conservative 
properties are suggested and demonstrated in the recent work [T3] . 

It seems rather difficult to construct numerical schemes for which one can prove 
the convergence to a (non-smooth) solution of the Camassa-Holm equation. This 
statement is particularly accurate in the case of general initial data and peakon- 
antipeakon interactions. Indeed, in this context we are only aware of the recent work 
[T2j in which we prove convergence of a tailored semi-discrete difference scheme to 
a dissipative weak solution. Before we can outline this scheme, let us discretize the 
spatial domain R by specifying the mesh points Xj = jAx, Xj^i/2 = (j + l/2)Ax, 
j — 0,±1,±2,..., where Aa; > is the length between two consecutive mesh 
points (the spatial discretization parameter). Let D, and denote the 

corresponding backward, central, and forward difference operators, respectively. 
The scheme proposed in [12], which is based on the formulation (II. 3|) . reads 

4-WJ + 1/2 + {uj + 1/2 V 0) + + {uj+1/2 A 0) + + D+Pj = 0, 

(1.7) ^ 

- D^D+Pj + P, = (uj+1/2 V 0)' + (uj_i/2 A 0)' + - {D-Uj+^/^f , 

where 

Uj+i/2it) ~ u{t, Xj+1/2), Pj{t) ~ P{t,Xj), for i > and j e Z. 
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If we interpret the Camassa-Holm equation ()1.4p as a "perturbation" of the 
inviscid Burgers equation, then the u-part of ()1.7|) might not come across as a 
reasonable (upwind) difference scheme. On the other hand, as pointed out in [12j . 
the key point is that with (jl.7p the quantity qj := D-Uj_^_i/2 satisfies a difference 
scheme which contains proper upwinding of the transport term in the equation for 
q := dxU, which reads dtq + udxQ + ^ + -P — = 0. Consequently, as is proved in 
[T^ , the scheme (|1.7p satisfies a total energy inequality, in which only the g-part of 
the total energy is dissipated (not the u-part, which is after all continuous). This 
is the essential starting point for the entire convergence analysis in |12| . 

The "semi-discrete" equation in ()1.7|) constitutes an infinite system of ordinary 
differential equations which must be solved by some numerical method. The main 
purpose of the present paper is to show that a fully discrete version of the scheme 
used in ^12J produces a convergent sequence of approximate solutions, and that 
the limit is a dissipative weak solution to (|f .ip . The fully discrete version that 
is analyzed in this paper is based on replacing the time derivative in (jl.7p by a 
forward difference, i.e., 

U, + l/2{t) ^ ^+"j + l/2 := , 

and evaluating the rest of (|I.7p at := nAt. Now u^j^i^2 should approximate the 
exact solution u at the point (i", Xj_|_i/2). This gives the fully discrete scheme 

^>;Vi/2 + (";Vi/2 V o) z?_u;Vi/2 + {u]+,f2 A o) i?+«;Vi/2 + = 

-D-D+P;^ + P; = (u^_,/2 V O)' + [u-^,/2 A O)' + i (^--"+1/2)' , 

where P" approximates P{t",Xj). As in [I2j this is a difference scheme which is 
tailored so that it gives an upwind scheme for the equation satisfied by q := dxU. 

The main aim of this paper is prove that the fully discrete (explicit) scheme 
(|I.8p converges to a dissipative weak solution of the Camassa-Holm equation. The 
starting point of the analysis is a total energy estimate, showing that the norm 
of the approximate solutions is (almost) nonincreasing in time. To this end, we 
must assume that 

(1.9) At^O{Ax^log{l + Ax^)) 

for some 9 > as Ax 0. This is a very severe condition, and it may seem that 
when using this method in practice one should use very small time steps. However, 
this is not a Neumann type stability criterion, and we do not have blow up if 
it is violated. Indeed, practical experiments indicate stability and convergence if 
At = 0{Ax). 

By appropriately extending the difference solution (jl.Sp to a function U/^xitj^) 
defined at all points {t, x) in the domain, we prove under condition (|1.9p that 
{''^Ax}ax>o converges strongly in to a dissipative weak solution of the Camassa- 
Holm equation p.ip - (|f .2p . Regarding the proof, we adapt the "renormalization" 
approach used in for the semi-discrete scheme, but there are several essential 
deviations and many parts of the convergence proof are substantially more involved 
and/or different. These differences are mainly due to the fact that the semi-discrete 
scheme, when viewed as a fully discrete scheme with "infinitely small time steps" , 
has a large and stabilizing numerical viscosity. Regarding the fully discrete (explict) 
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scheme (|1.8|) , to account for this lack of numerical viscosity the convergence analysis 
relies heavily on the CFL condition ()1.9p and differennt With reference to the 
differences between the semi-discrete and fully discrete schemes, let us here point 
out just one aspect, namely that the norm of the fully discrete approximation is 
not entirely nonincreasing but can grow slightly with a growth factor that, however, 
tends to zero as Ax 0. Compared to semi-implcit case [T2l, the proof is notably 
more complicated and involves working with a version of the scheme (|1.8|) in which 
the quadratic terms have been suitably truncated. 

The paper is organized as follows: In Section [5] we introduce some notation 
to improve the readability and recall a few mathematical results relevant for the 
convergence analysis. The finite difference scheme and its convergence theorem are 
stated in Section [3l The convergence theorem is a consequence of the results proved 
in Sections |4][8l Finally, we present a numerical example in Section [9l 

Throughout this paper we use C to denote a generic constant; the actual value 
of C may change from one line to the next in a calcuation. We also use the 
notation that < bi to mean that < Cbi for some positive constant C which is 
independent of i. 



2. Preliminaries 



In what follows. Ax and At denote two small positive numbers. Unless otherwise 
stated, the indices j and n will run over Z and 0, . . . A^, respectively, where NAt = T 
for a fixed final time T > 0. For such indices we set Xj = jAx, Xjj^ij^ = (j + l/2)Ax, 

= nAi, and introduce the grid cells 

= [x,_i/2,a;,+i/2), = and /; - x /". 

The following notations will be used frequently: 



a V = max {a, 0} 



a + \a\ 



a A = min {a, 0} 



a — \a\ 



For n e {0, ...,Ar}, let = denote an arbitrary sequence, where n 

refers to "time" and j to "space". We will frequently employ the following finite 
difference operators: 



:= 



Ax 



D,v- := 



3-1 



Ax 



J • 2 
We also use the notations 



2Ax 



and D]_v!l 



At 



I < p < oo. 



sup I Vj 
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Occasionally, we also use the "space-time" norms of v = {^'"In^o ~ {^i }jn' 

i/p 



\ n=0 I 



Note that \i v ^ P , p < oo, then limj^±oo Vj = 0. 

Let {vj}j^2, ^'^d {wj}j^^ denote two arbitrary (spatial) sequences. Suppose 

< oo. Then the following discrete Sobolev inequality holds: 



hi 



(2.1) 



< 



1 

V2 



The discrete product rule takes the form 

(2.2) D± (vjWj) — VjD±Wj + D±VjWj±i. 

Moreover, the discrete chain rule states 



(2.3) D^fiv,) = nv,)D±v, ± ^f"{^f){D^v,)\ 



fee' 



for some number between Vj±i and Vj. 

We continue to collect some handy results for later use, starting with a discrete 
Gronwall inequality. 

Lemma 2.1. Assume that c'^ > and > Q for all k — 0, . . . , N , and that the 
sequence {w"}^^o satisfies the difference inequality 

(2.4) Z)*^?/" + /" < c"w", n = 0, ...,7V-1. 

// u" > for alln = 0,...,N, then 



N-l 



N-1 



N-1 



exp At J2 c" At J2 - ^''P H ^" 



n=0 



r!=0 



fc=0 



ri=0 



Proof. Set i?" = cxp (^-At X^Lo j ■ Then we have 

1 



At 
exp 



(-Atf]cMi^(l-exp(Ate")) 

< (1 - (1 + c"At)) = ~c"i?"+i. 

Hence, multiplying (|2.4p with we arrive at 

Multiplying this by At and summing over n, we see that the lemma holds. 



□ 



The next lemma contains estimates for the solution of a discrete version of the 
differential equation P — d^^P — f. 
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Lemma 2.2. Let {fj}j^i he a sequence in fl i'^, and denote by {Pj}j^2, 
solution to the difference equation 

(2.5) P,-D^D+P,^f,, jeZ. 

Introducing the notations 



h= 1 + 2 



1 - e- 



, , , Ax Ax r — J 

K = In ( 1 + + — \/4 + Aa;^ 



(Aa;): 

the solution {Pj}j^^ takes the form 

(2.6) P,. =/i^e-«l^-V., jeZ. 



Moreover, the following estimates hold: 
(2.7) 

(2.8) 



(2.9) 



{PA, 



< C 

< C 



where C > is a constant independent of Ax. 

Proof. To verify the solution formula (|2.6p . we define pi by 

for some constants c and cr yet to be found. We shall choose these so that 

'l, ifi = 0, 



0, otherwise. 



{I-D+D-)p,= 



If we find that this holds with a ~ k and c ~ h then (|2.6p holds. We observe that 
for i ^ 

I*locosh(a)-1 



D+D-p^ = ce-'"i*i2 



Aa;2 



Hence cr must satisfy 



a — cosh ^ I 1 

For i = we find that 

Po - D+D-Po c I 1 



Ax^ 
2 

A? 



-1) 



If this is to be equal 1 then c — h. 

For later use, one should observe that 

(2.10) h=^+0{Ax% ^^ = l + 0(Ax), \£^-li = l + 0{Ax). 
For any j £ Z, we have 

Furthermore, 



i/d<iimii,i 
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Hence, we have proved (|2.7p 
From ([2^1) . 



Ax 

g-K|i-j-l| _ g-K|j-i| 



Aa; 



hT.- /^ + ^ E ^ — 



Aa; 



i— J -^— — 00 

Using (|2.1Qp we acquire from this the following two estimates: 
|Z?+P,|</i^e-«l-^-||/,|<||{/0||,. 



and 



{D+P,}^ < hAxJ2e--\^-^\ |/,| < ||{/.}||,, , 



Therefore holds. 

It remains to prove (|2.9[) . To this end, we multiply the equation (|2.5[) by AxPj 
and perform a summation by parts to discover 

2 1 



from which (12.91) follows 



□ 



We shall routinely use some well-known results related to weak convergence, 
which we collect in a lemma (for proofs, see, e.g., [24]). Throughout the paper we 
use overbars to denote weak limits. 

Lemma 2.3. Let O be a bounded open subset o/R*^, with M > 1. 
Let {i'n}„>]^ be a sequence of measurable functions on O for which 

sup / <^{\vn{y)\) dy < oo, 

n>l Jo 

for some given continuous function $ : [0, oo) — > [0, oo). Then along a subsequence 
as n —> oo 

g{vn) g{v) in L^{0) 
for all continuous functions g : M — )■ R satisfying 

l^l^oo $(|w|) 

Let f;: R — > (— oo,cx)] be a lower semicontinuous convex function and {wn},j>]^ a 
sequence of measurable functions on O, for which 

Vn ^ V in L^{0), g{vn) G L^{0) for each n, g{vn) g{v) in L^(0). 

Then 



liv) < gi''^) fl-e- on O. 
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Moreover, g{v) G L^iO) and 

/ g{v)dy <\immi / g{vn)dy. 
Jo Jo 

If, in addition, g is strictly convex on an open interval (a, b) C 



g{v) — g{v) a.e. on O, 

then, passing to a subsequence if necessary, 

Vn{y) v{y) for a.e. j/ G {y e O | v{y) G (a, 6)}. 

Let X be a Banach space and denote by X* its dual. The space X* equipped 
with the weak-* topology is denoted by -^^^^cak' while X equipped with the weak 
topology is denoted by Xweak- By the Banach- Alaoglu theorem, a bounded ball in 
X* is (7(X*, X)-compact. If X separable, then the weak-Tk- topology is metrizable 
on bounded sets in X* , and thus one can consider the metric space C ([0, T]; X^^^j^) 
of functions v : [Q ,T] ^ X* that are continuous with respect to the weak topology. 
We have ^ in C {[Q ,T\] X^^^^) if {vn{t) , (t>) x\x {v{t),(j))x*,x uniformly 
with respect to t, for any <f) € X. The following lemma is a consequence of the 
Arzela-Ascoli theorem: 

Lemma 2.4. Let X be a separable Banach space, and suppose w„ : [0,T] X* , 
n = 1,2, is a sequence of measurable functions such that 

ll"n|lL°=([0,T];X*) - ^' 

for some constant C independent of n. Suppose the sequence 
[0,T]3t^ {v„{t),^)x*,x, n = l,2,..., 

is equi- continuous for every $ that belongs to a dense .subset of X. Then w„ 
belongs to C {[0,T]; X^^^y.) for every n = 1,2,..., and there exists a function 
V € C {[0,T]; X^^^^) such that along a subsequence as n —> oo 

v,,^v mC{[0,T];Xt,^^). 

3. Explicit scheme and main result 

In this section we present the fully discrete (explicit) difference scheme for the 
Cammassa-Holm equation (|1.3p . which generates sequences |u"+i/2| and {-P"} 

for {n,j) G {0, . . . , A^} x Z. We let |it"_^i/2} solve the explicit difference equation 

(3.1) i?V""+i/2 + (""+1/2 V O) + {u]+,/2 A O) i?+^^;Vl/2 
where the initial values are specified as follows: 

(3.2) U°^i/2 = "o(-i^j+i/2), 
Given |w"^j^^2|! determine {-P"} by solving 

(3.3) - D^D+P;^ + P; = (^.;Vl/2 V O)' + (u]_,^, ^ + ^ {0-^^1/2)' 



D+P^ = 0, 



which is a linear system of equations that can be solved as outlined in Lemma l2.2l 
Next, let us derive the difference scheme satisfied by 

(3.4) = i?-w;'+i/2- 
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This will be done by applying the difference operator D- to the u-equation (j3.ip 
To this end, we apply the discrete product rule to find 



and 

so that 
(3.5) 



J + l/2 



J + 1/2 



'■3 + 1/2 



A £> I u 



+ "j + l/2 



= (""-1/2 V O) D^q, + (u-^,,^ V O) 
= (""+1/2 A O) D+q, + («;Vi/2 A O) g,, 



.;Vi/2 V O) D^u-^„, + [u]^,/, A O) D+u]^,/,_ 

= (""-1/2 V O) D^q^ + (u]^,^, A O) + (?,")' • 



The P-equation (j3.3[) rephrased in terms of g reads 

2 



- D^D+P" 



^j + l/2 



vo 



'J-1/2 



AO 



2^. 



(3.6) 

Employing (|3.5p and (|3.6p when applying Z3_ to the u-equation in (|3.ip yields 
I^V-?; + (""-1/2 V O) i?_g; + (ii;Vi/2 A o) 



(3.7) 



pn 



"5Vi/2 V 



""-1/2 A 



0. 



2 ' 'J 

Regarding the initial values, in view of (|3.4p and (|3.2p . we observe that 



(3.8) 



Since the variable q = d^u can be discontinuous, (|3.7p represents a natural 



upwind discretization of the equation for q, dtq + udxq - 



P = 0. 



The main result of this paper is the convergence of the scheme to a dissipative 
weak solution of (|l.ip - p.2p . which is defined in the following sense [39| 1401 : 

Definition 3.1. Fix a final time T > 0. We call a function u: [0,T] x R ^ R a 
weak solution of the Cauchy problem for (|l.ip - (jl.2p on [0,T] x R if 
(D.l) MeC([0,T] xR)nL°°(0,r;7Ji(R)); 

(D.2) For all s and t in [0,T], with s <t, we have \\u{t, •)IIhi(r) ^ ll"('5; ')lli/i(R)- 

(D-3) u satisfies (|1.3p in the sense of distributions on (0,r) x R; 
(D.4) u{0,x) ~ uq{x) for every x e R; 

(D.5) If in addition, there exists a positive constant K such that 



2 

Ux{t,x) < - + K ||uo|Ihi( 



{t,x) e (o,r] X R, 

then we call u a dissipative weak solution of the Cauchy problem (|l.ip - (jl.2p . 



In addition to dxU E i°°(0, T; i^(R)), cf. (D.2) the dissipative weak solutions 
u that we construct in this paper will possess an improved integrability property, 
namely d^u S Lf„^((0,T) x R) for p < 3, i.e.. 



\dxuf dxdt < C{a,b,T,p), Va,6 G R, a < b. 
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To state our main convergence result and also for later use, we need to introduce 
some functions (interpolations of the difference approximations) that are defined at 
all points (f, x) in the domain. We begin by defining the functions 



u,+i/2{t) = 7i;'+i/2 + {t- n D\u'hi/-2, 
With the aid of these we define 

(3.9) qc,{t,x)=qj{t), {t,x)&I^ 
and 

(3.10) UAx{t,x) = Uj^i/2{t) + (x - Xj_i/2) qj{t), for (i, x) e 

for j S Z, n = 0, . . . , TV — 1. Note that 1 UAx{t, x) is a continuous function, since 
Uj_i/2(t) and qj(t) are continuous. Regarding the continuity in x we have that 

lim UAx(t,x) ^ Uj-i/2(i) + [xx^iji - Xj-1/2) Qjit) 



for t e /" 



- ^-1/2 + it- n IA^;-U2 + Ax {q- + {t- n D\q-) 

= ""-1/2 + n ^v<-i/2 + (";Vi/2 - ""-1/2) 

+ it- n D\ (u;Vi/2 - ""-1/2) = ",+1/2 (t), 



and therefore uax is continuous, and furthermore dxUAx — qAx almost everywhere. 
Observe also that, due to (|3.8|) . there holds qAx{0,x) dxUo in L^(R) as Aa; 0. 
Similarly to uax, we define a function Pax by bilinear interpolation. First, let 

p^.(t) = pj' + {t- t")D'^P^, t e 

and then define 

(3.11) P^.x{t,x) = Pj{t) + {x~ Xj)D+Pj{t), {t,x) e 

for j e Z, n = 0,...,iV- 1. 

We are now in a position to state our main result. 

Theorem 3.1. Suppose (|1.2p holds. Let {uax}ax>o sequence defined by p.lOp 
and (|3.ip - p.4p . Then, along a subsequence as Ax J, 0, 

UAx^u ini/i„,((0,r) xK), 

where u is a dissipative weak solution of the Cauchy problem (jl.ip - (jl.2p . 

This theorem is a consequence of the results stated and proved in Sections |3][5] 

4. Total energy estimate and some consequences 

The purpose of this section is establish a discrete total energy estimate for the 
difference scheme (I3.1I)-(I3.4I). 



Lemma 4.1. Assume that Ax and At are related through the CFL type condition 

log(l + A/) Ax^ 



(4.1) At < 



12 



G. M. COCLITE, K. H. KARLSEN, AND N. H. RISEBRO 



where C is a constant (to he detailed in the proof of the lemma) that is independent 
of Ax and uq and 6 > 0. Then, for any Nq G {0, . . . , N}, and for all sufficiently 
small Ax, 



(4.2) \\u^o\\^^,+Ax'AtY^Y^ 



^j + l/2 



Proof. For the proof of (|4.2p . we shall need to introduce an auxiliary difference 
scheme. To this end, we start by defining the cut-off function 

{-M, u<M, 
ue[-M,M], 
+M u > M, 

where A4 > is a fixed constant (to be determined later on). Now, let |'ti"_^i/2| 
and l-P/'l solve the following system of difference equations: 



D 



(4.3) 



+^;vi/2 + (/''(^;'+i/2) vo) D^u-^,/, 

for n = 0, . . . , iV — f and j £ Z, and 

(";Vi/2) V o) (f.;Vi/2 V o) 

+ {f" (-"-1/2) A o) (s;U/2A0 



V(r^(^^^,/2))^-""+i/2, 



2 

for n — 0, . . . ,N and j e Z. Regarding the initial data, we set — w° for j S Z. 

If we define g" D^u^^^^^, then it is straightforward to see that satisfies 
the difference equation 



(4.4) 



U7 + (-"-1/2) V 0) D^q^ + (/^^ A 0) D+q_ 

+ 1^- (/" (-"+1/2))^; 

- {.r («;vv2) V 0) V 0) 

- (""^1/2) A 0) {u-_,/,A0)+Pp=Q. 



Multiplying (|4.3p by w"_^_^^2 '^^ f^^^ that 

«^+l/2i5^""+l/2 + (-;Vl/2) V 0) V 0) 

+ (/ [^1+1/2) A 0) (^^+1/2 A 0) g^i + ^+i^"s;Vi/2 = 0, 
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while multiplying ()4.4p by g" gives us 

q-D\q- + (^^--1/2) V 0) {D^q-)q- + [f'' (u;^!/^) A o) {D+q-)q 



(4.6) 



(-;Vi/2) V 0) [u^^y, V 0) 



/^(-^i/2)A0 



A0)g"+P"q7 = 0. 



Adding (|4.5p and (|4.6p . multiplying the result with Aa;, and summing over j yields 
Ax 5^ (z2j+i/2i?^uJ+i/2 + q]D\q^) + I + II + III = 0, 

where 

+ A. E (/" A 0) + ^ E ^- (/ 

II = AxE (^"+1/2) V 0) V 0) 

+ Ax E (-;Vi/2) A 0) (u-^,i, A 0) g;v, 

iez 

-AxE(r (-;W2)vo) (^;w2vo)g; 



"'S + l/2))(9") 



Ax E ('^i-1/2) A 0) (u"_i/2 A 0) = (by shifting indices), 



III = Ax E ^+^"""+1/2 + Ax E -P"?" = (by summation by parts, of. 

j6Z jez 



Let us now deal with term I. The discrete chain rule 



tells us that 



Ax. 
2" 



T^{D^~q-f. 



Hence 



I=AxE(/''(^--i/2)vo) 

+ Ax E -+1/2) A 0) 



iez 



Ax 



(97)' = Ii+l2, 
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where 



Ii = Ax^(r^ (""-1/2) VO) 



Ax' 



■E[(/"' (-;Vi/2) VO) (i^+g}')^ - (/^^ (^;Vi/2) Ao) 

■E|/"fei/2) (^+^")'>o. 



To handle the Ii-term, we use the discrete product rule 





L'4 



(r'(^"-i/2)vo 
(/" (""-1/2) A 
= {f' (^;vi/2) A 
= {f' (";vi/2) A 

Using this we find that 

Ii = A^E^- 



'j+1/2 



vo 



r^(-"-i/2)Ao) 

f ' (""+1/2) A 0) 



'J + 1/2 



vo 



AxJ2D-(f 
Ax^D 



^j+1/2 



VO 



(r^(""-i/2)A0) 



AxE^- (/*^ (^7+1/2) A 0) 
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-A.E^-(/"(""+i/0) 



Ax 



0. 



Summarizing our findings so far: 

7"ez 

(4.7) 



J + 1/2 



^+9- =0- 



Next, by (f^ . 



(4.8) 



Ax 



E ((-^Vi/2)^ + ('TjT) 

iAtAxEfteVi/2)' + (i?U;) 



Hence, we must now estimate 



AiAxE(K-;Vi/2)' + (i?v9;)') 



Using (|4.3p . (14. 4p . and the basic inequality ^X^fci*^') — 2' ^ X]f=i('^^)^' which 

holds for any sequence {a^}^^]^ of positive real numbers, there is a positive constant 
ci that does not depend on Ax such that 



AiAxE ( 

< ciAiA^E 



(r(^;W2)vo)^(^7)^ 



7+1/2 AO (g^Vi) + 



'j-l/2 



VO) (D-g 

2 



/"(""+i/2)))'(€) 



'i+i/2 



vol (w;Vi/2V0 
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+ (/^ (^"-1/2) A 0)' (u]_,^, Ao)\ (p^y 

< CiAt {Ji + J2 + J3) , 



where 



Ji = Aa; ^ 



J2 = Ax ^ 



J3 = Ax^ 

fM 



{r {u]^,/.) V 0)' {q^y + {u]^,/.) A 0)' (g^Vi) 

+ (/^ («^+i/2)vo)'(n^^,/,vo)' 

+ (/^ A 0)'(^^_,/, A 0) 

(/"(-i-i/2)v0)'(l?_C)' 
+ (/- (.,+v2) A 0)^ (D+^»)^ + {d. (r 



Since / (u) < M, the following bounds hold: 



(4.9) J,<2M'AxY^ +{Q7f 



and 



(4.10) 



(/''(^"-i/2)vo)'(2?_g;Ax)' 

+ (/*'(«"+i/2) A0)'(i?+g7A.^^ 

(/-(«^+,/,))a.)'(,7)^ 



+ 



■Aar^ 



for some constant C2 > independent of Aa;. 
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To estimate J3 we use Lemma [2.21 specifically p.9p . which implies that 



J3 < csAx^ 



VO) (m^+i/2V0 



(4.11) 



(/"(""+l/2)))"('^;f 



< C3M2 I 1 



for some constant C3 > independent of Aa;. 

Blending (|49)) . (|4.10p . and (|4.1ip we derive the bound 



At 
A? 



(4.12) AtAxY, (^V«;'+i/2) + (^V^^O' < CM' [At 

where the constant C is independent of Ax. 

Combining (g^l), and ((i7T^ . it follows that |u"+i/2} obeys the following 

discrete energy estimate: 



t ll,--,n||2 



Ax' 



+ II" llhl 



< CM' ( A< 



At 
A? 



By the discrete Gronwall inequality, cf. Lemma [2.11 



No~l 
n=0 



Choosing M = ||ito||jLfi(-jj-, and recalling the CFL type condition (14. ip . we deduce 



< (1 + Aa;") 



e\NAt 



< e 



t" Ax" 



< 2 if Ax*' < 



log 2 
T ■ 



Therefore, in particular 

||i2"||^i < V2||?/°||^, forn = 0,...,iVo. 
By the discrete Sobolev inquality (|2.ip . we find that 

< i= \\u^^, < < M forn - 0, . . . , / 

This means that u" will "never notice" /*^, since 

(""+1/2) = ""+1/2 for j e Z and n = 0, . . . , iVo 
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Therefore, 



u]+i/2 = """+1/2 and = for j e Z and n = 0, . . . , TVo- 



Finally, (|4.2I) follows by noting that 



gCM^^(t"o-i_r) > ^ for n = 0, . . . , TVo - 1. 



We conclude this section by stating some immediate consequences of (|4.2 
Lemma 4.2. For n = 0, . . . , N , 



□ 



m^,n^<c\\uo\\l^^ 



\\D+P^,^,\\D+P^,^<C\\uo\\H^^^), 

where C > is a constant independent of Ax. 

Proof. This follows immediately from Lemma |2.2[ noting that in this case 

= (<+i/2 V o)' + ^^y + l 



and thus 



|/,IU<h"||^<(l + Aa;)^ 



\u 



□ 



5. One-sided sup-norm estimate 



Lemma 5.1. Assume that At satisfies the CFL type condition (j4.ip and that Aa; 
is sufficiently small. For n — 0, . . . , N and j G Z, we then have 

(5-1) <i^<§: + ch'\L^ 

where C > is a finite constant. 

Proof. We can write the difference equation for {q"}, see p.7p . as 

q^+' - q'l (1 - Aa - A6) + q]_,Xa + q^;^,\h - At^ 

1 



+ At(^(w;_i/2V0)' + (uJ+i/2A0 



pn 



where 



a = A (m"_i/2 V , 5 = -A (ujVi/2 ^ o) . = At/Ax. 
Now we have uniform bounds on 1 1 u" I Loo and and thus 



(5.2) g^+i < (1 - Aa - \h) + gJ^iAa + q'^^^Xh - At^ + AtL, 



for some finite constant L < llu'^llf , 



Set g" = max {<??, We claim that 



2 



19. 
2 

if Ai is chosen sufficiently small. 



(5.3) g^+i < q'; - Ai^^ -f AtL, 



AN EXPLICIT SCHEME FOR THE CAMASSA-HOLM EQUATION 



19 



First we choose At so small that A(a + b) < 1/2. Then if qj = qj the claim 
follows immediately from (j5.2|) . 

Next, assume = q^j^i- Then note that 

('^;)' = ('?;-i)' + ('^;)'-('^^i)' 

^{q^y + '-^^A^D^q^ 

= {,Ur + "-^^^^^^^D^q-. 
Since D^q^ < 0, we find that 

Using this we can rephrase (|5.2p as 

< 9" (1 - A (« + + b)) + q^;_,X (a + + q^^.Xb 

_ At^^^ + AtL 

< q^ - At^ + AtL F (q^) , 

if A < i. The proof of ([531) if = is similar. 

Note that F'(q) = 1 — Atq, and thus F is increasing for q < 1/ At. Furthermore, 
by the CFL type condition (|iTT|l . At = 0{Ax^), and by the bounds on ||M"i|^cx>, 



\AxJ \At^/V ~ At' 
for sufficiently small At. Therefore, setting A/" — maxj g", from (|5.3p we get 

M"+i < F (M") . 
Now set = M" - \/2l. Then 

^n+i < f + ^211 - V2L 



= + %/2Z- ^ (z" + \/2Z)^ + LAt - V2L 
= (l - AtV2Lj - ^ {Z'^f 

Now, clearly if Z" < 0, then < 0. Hence, if Z° < 0, then Z" < for aU 

n > 0. If Z" > 0, by [Ml page 271], 

2 2 

Z" < < _ 

- + 1/ZO - i"' 

This finishes the proof. □ 



^This maximum exists since G ■ 
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6. Higher integrability estimate 

We begin this section by deriving a "renormalized form" of the finite difference 
scheme for ^j, so let / be a nonhnear function (renormahzation) of appropriate 
regularity and growth. Multiplying (|3.7p by f'{qj) and using the discrete chain 
rule, which in the present context reads 

where ^"±1/2 ^ number between qj and <?j±i, and gj^^^^ is a number between 
and q"^^. Multiplying the scheme (13. 7|) with f'{q^) we obtain 



Dif {q-) + [u]_,/, V 0) D^fiq-) + A o) D+f{q^) 



(6.1) 



pri 



/" ( 



// / n+1/2 



where 



IaxJ",] := 



Ax 



-("■+1/2A0) nqJ^,/,)iD+q-y 



Let us now write ()6.ip in divergence-form. To this end, observe that the discrete 
product rule (|2.2p implies the following relations: 



(";vi/2vo)/(g;) 

"-i/2A0)/(g7) 



V 0) D^fiq^) + (u^^,/^ V 0) fiq^), 
;Vi/2 A 0) D+f{q^) + D+ (u'J_,/^ A 0) fiq^) 
(u;Vi/2 A 0) D+fiq!^) + («;Vi/2 A 0) /(g;). 



and therefore, using that g" = ^-"j 



V 0) i?-/(g;') + (m;Vi/2 a 0) D+fiq^D 



Vi/2 V 0) fiq^J)] + D+ [(z.;Li/2 a 0) fiq^)] - g;/(g,"). 
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Hence, we end up with the foUowing divergence-form variant of the renormahzed 
difference scheme (16.11): 



D\f{qY) + V O) f{qj)\ + D+ [{u-^^/2 A o) f{q]) 



(6.2) 



Q 



// / n+1/2 



To ensure that the limit of (OxUax)^, cf. (|4.2p . is nonsingular (i.e., not a measure), 
we shah need the following higher integrability estimate: 



Lemma 6.1. Let q/^x be defined by (|3.9p . and assume that the CFL type condition 
()4.ip holds. Then, for all finite numbers a, b, a with a < b and a € (0, 9), 

i-T rb 



(6.3) 



/ / \qAx 

Jo J a 



2+a 



dx dt < C, 



for some constant C = C{a, b, T, a) that is independent of Ax. 



Proof. Define 77e(g) = y^e^ + q^ — e. Note that rj^{q) \q\ for small e and 



-1 < v'Aq) 



< 1 and < 77" (q) 



(£2+^2)3/2 



< 



V a/2 



In (|6.2p . we then specify 

f{q)^Ve{q){l + q 
One can easily check that 

fiq) = vUq) (1 + q'f" + c^vMq (i + qT^'^' 

fiq) = <(9) (1 + qT^' + 2ar^'Mq (l + qT^'^' 

+ ar,,iq)^{l + qY'-'+2{^-l) {l + q') 
SO that in particular f"{q) > and 



q/2-2 



(6.4) 

Next set 



nq)^rj':iq){l + q') 



2\"/2 



bounded terms. 



H{q) ■■=\.f\q)-qf{q) 

= ^ (v'M) (1 + qT" + -VM (1 + qT^"'' q) qvM (l 

/ wi , 2^"/2 \qveiq) 
= qveiq) (1 + g 

= We (<?)(! + 9' 



2\"/2 



2^a/2 1 



277,(g) 2(1 + ^2) 
aq-^ e 



- 1 



- 1 



= : HMheiq), 
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with 

Now note that 



a/2 



he{q) = 2 



aq 



- 1 



a - 1 

hm hi,{q) < lim hi{q) = — - — < 0. 

q — ^ — 00 q — > — 00 2 

Hence, for a,e < 1, we can find a constant K > such that 
(6.5) 



h^{q) < for aU q < -K. 



Let us continue by defining the sets 

AA" = {j e Z I < ~K} , C^ = {3eZ \ -K< q^ < 0} , 

and V = {j el I > 0} , 

where K is defined in (|6.5p . Moreover, let < x{x) < 1 be a smooth cutoff function 
satisfying 

{0, a; < a — 1, 
1, xe [a, 6], 
0, x>h+l. 

We multiply (|6.2[) by x(a;j)AiAa; and sum over (n, j) € {0, . . . , — f } x Z to get 



(6.6) 
(6.7) 

(6.8) 

(6.9) 
(6.10) 

(6.11) 



< -AtAa; J] J] (g^) x (x,) 

- AtAx^ ^ /i, (g«) ff, (g^) X (x,) 
-AtAxY,[D^ ((u]^,/, V 0) / (<z;)) x(2;,) 

^+((^;-i/2A0)/(g7))x(x,) 



AtAa;^A;'x(a:,) 

x(x,) {D*:,q-y 



where 
(6.12) 

Now for j E A/"" we have that 



SVi/2 V 



^^1/2 A 



1 — a 



1/2 



< 
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Therefore 

< |(I611)| + |(I6I1)| + |(I611)| + |(I61])| + |(I^ 

We shall now find bounds on all the terms on the right hand side; in what follows, 
we let C denote a generic constant independent of Ax, s, and a. 
We start with ([OT|) . By ([XT]) 



(0 



VO D-g' 



^j + l/2 



A 



+ (pry + (<+i/2 V o) V {u]_,/, A o) V (,;)^ 



To bound the "integrals" of these terms we must use the CFL type condition (|4.ip . 
which imphes that A< = 0(Ax^+''). First, 

i 

< CAx'+' J2 i^^D^q^f < CAx'+'> {q]Y < CAx'. 

3 3 

Therefore 

AtAa;^ At ((<_i/2 V o) D^q"])" + ((h;Vi/2 A o) D+q^^y x{x, 

n,3 L 

< CAx'^T -> as Ax -> 0. 
We also find that 



axy: 



3 L 



f + (";'+i/2vo)V(^i;'_i/2^o)^ 



< CAxE 

3 

< CAxE 



LP" 



LP" 



+ (";'+i/2V0)'+(zi;Li/2A0 

+ (";'+i/2) 



since uax and P" are uniformly bounded. Thus 

AtAxY^At \{p-j' + («;Vi/2 V O) V iu^_,/, A O)' 



as Aa; 0. Additionally, 

Ao; E At (9;')' < Y (A:e9,")' (^^O 
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<CAxi+''^(q7)'<CAa;^ 



AxAt^ At {q"Yxixj) < CTAx^ -> as Aa; -> 0. 



Now we have estabhshed that 

(6.13) AxAt^At{D'i^q^yx{x3) = o(^Ax'''^ as Aa 



X -> 0. 



Recalhng (|6.4p this imphes that 
IdSUDi < CAxAtY, — fl + U 



n+l/2 



a/2 



{D\q-)\{x,)+o[Ax') 



When we estabhshed (|6.13p we always had a "Ax^ to spare" , which we can use 
now. With /3 = — a > 0, we get 

< CAxAt^ ^Ax" (l + {qj^'^^yy Ax' {D\q^;f 

< CAxAtY, ^ (ax^ + (Axg;+^/')T^' Ax' {D\q^Y 



< C^AxAtY^Ax' {D\q^f < C^, 



since 



AxV (Axg;+^/^)y^' < (Ax^ + 4 



(Ik" 



a/2 



< c. 



Now choosing e ~ Ax^ ^ we finally conclude that | (|6.1ip | is bounded. 

Next we turn to (16. 6p . For —K < q < 0, we have that \h^{q)H^{q)\ < C, where 
C is independent of e. Therefore 

\&\ < CT{b~a + 2). 

To estimate | (16. 7p | , observe that 

\m\<cAtAxj2 E + 

Now, equipped with (|4.2p and (|5.ip . it is possible to bound the right-hand side by 
a "Ax independent" constant exactly as was done in [12] . 
Regarding (|6.8p . observe that 



AtAx E ^- [(<+i/2 V O) / (g^)] x(jAx) 

AiAxE (""+1/2 V O) D+xijAx)f {q^) 
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<Ct|^sup||w"||,^ \{D+x{xj)}^ 
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+ sup 



{k"|-}, 



{D+xii^x)}.. 



cT |^sup||u" 



{D+x{xj)}^ 



supljg 



{D+xU^x)}.^ 



Therefore, also | (|6.8p | is bounded independently of Ax. 

Next we focus on (I6.9p . Remembering that |/'(q)| < C(l + we find 

»| <AtAa;^(||P"||,^,+2h"||,2^) \f'iq'^)\x{xj) 

< CAtAxJ2 {WP'le^ + 2 (1 + \q^\") x{xj) 

2 sup 1 1 u" 

n 

+ sup 



< CT sup IIP" 



<CT sup||P"||^^+2sup||u"||;„ 



^1 



+ sup||g" 



{xO'Aa;)}, 



{x(jAx)}^. 

Finally, keeping in mind that / > 0, we treat (|6.10p as follows 

|MI<Ax^(/ (gf)+/(9°))x(^,) 
j 

<CAxJ2(2+\qff + \q^\-)x{x,) 



< 



I iV||l+" I II 0||1+" 



Summarizing, we have established 



1+ 

{xixj)}j 



a/2 



AtAxj2Y.h-\v.{q]){i + {q-y) <c. 

n j 

The statement of the lemma follows by noting that 

\qr-<\q\Uq){l + qr'' + \q\'^\ 
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and using, in combination with (|4.2p . the bound 



□ 



7. Basic convergence results 



The purpose of this section is to present some straightforward consequences of 
the a priori estimates estabhshed in the foregoing sections. More precisely, we 
prove that the two sequences {uAx}/\^xyQ, cf. (|3.10|) . and {-Pax}ax>0' I|3-11P - 
have strongly converging subsequences, starting with the former. 

Lemma 7.1. There exists a limit function 

u(zL°° (0, T; H^{R)) n C([0, T] x R), 

such that along a subsequence as Ax — > 

(7.1) UAx^u tnL^{0,T;H\R)), 

(7.2) uax — > u uniformly in [a, b] x [0, T], for any a < b. 
Additionally, 

(7.3) t 1-^ ')IIhi(r) non-increasing, and 

(7.4) Vim u{t, x) ^ uo{x), xeR. 

Proof. First we note that for < G 

/ {uAx{t,x))^ dx^y2 -l—{Xj + i/2 - x)Uj_i/2{t) + {x - Xj_i/2)Uj + i/2{t) 



'2 



< ^ (^(* - Ax E (-;:i/2) ' + it"^' - Ax E (""+1/2) ' j 

md 

y" idxUAxit,x) f dx = [qjityf 

< ^ (|(t-nA.;E«+^)%(r+i-OAxE('?")'j • 



Hence, 



\\uAxit,-)\\l^M) < ^ {{t-n + (^"^' -0 



At 

Let s G 7™ with m < n. Then, using 



WuAxit, OIIh^W < ^ ((i - + (^"^' - ^) ll""ll?^MK)) 
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At 



+ (t"+l-<)(||u™||^^,(„)-hA.(.,-)ll?,.(M)) 



,m+l||2 



Lf/i(E) 



Lh"1(R) 



2 

This implies ([771]) and (TTSl) . 

Next we prove that {ua2;}ax>o uniformly bounded in W^'^+"((0, T) x (a, 5)). 
We can assume that a — 1 — Xj^ and b + I = Xj^ for some integers ja and jb. Since 
5 1-^ |g| IS convex, 



(7.5) 



for some constant C — C{a, a, b, uq), where we have also used ()6.3p 
Now set a — {x — Xj^i/2)/Ax. Then, for x £ Ij, we have 

\dtUAx{t,x)\ = (1 - o-)i:>+uJ_^/2 +^^+^"+1/2 

< (1 - (7 



^>"-l/2 



^+"i+i/2 



Furthermore, by the uniform bounds on D-^-P^^ and max, 



-^+"i+i/2 



<C(l + k"| + |g" I). 



Using this. 



- 7 

nb / N jt 

\dtUAx\^^°' dx <C \1 + AtJ2^xY^ \q^ 



2+a 



< c. 



Now {uax}ax>o C cc C"^^ on (0, T) x (a, 6) with £ 1 - 2/(2 + a). 

Therefore, along a subsequence, uax ^ uniformly in (0,T) x (a, 5) as Ax —^ 0. 

Let us show that the limit satisfies the initial condition (|7.4p . Fix a; S K and let 
t £ (0, 1). We have x G Ij for some j and UAx{xj^i/2,0) = uo(xj_i/2), so that 

\uAxiO,x) - Uoix)\ < \uAx{0,x) ~ UAx{0,Xj-i/2)\ + |uo(a^j-l/2) " Uoix)\ 
<C{x- Xj_i/2Y ■ 

Consequently, 

\u{t,x) - Uo{x)\ < \u{t,x) - UAxit,x)\ 

+ \uAx{t, X) - UAx{0, X)\ + \uAx{0, X) - Uo{x)\ 

< \u{t,x) - UAx{t,x) \ + Ct^ + CAx^. 
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Now we can let Ax and then t ^ to conclude that u{x,0) = uo{x). This 
draws to a close the proof of the lemma. □ 

Lemma 7.2. There exists a limit function 

P e L°°(0, T; W^'°^iR)) n L°°(0,T; W^-\R)), 
such that along a subsequence as Ax — > 

(7.6) Pa. ^P m LL((0,T) xR), 1 < p < cjo. 

Proof. By the bounds on P" in Lemma 14. 2[ we see that Pax is bounded in L°° 
uniformly in Ax. Next we show that {dtPAx}Ax>o bounded in L^{{0,T) x R). 
For t e [f\t"+^) and x G Ij+1/2, 

dtPAxit, x) = D\P'^ + {x~ Xj)D\D+P!^ = {l-<y) D\Pf + <jDXp;\i, 

where a = {x - Xj)/Ax. Write £)*^P" = X" + YJ\ where X" and Y," solve 



rr - D^D+Vr = -D 



-1/2 
2 



AO 



""+1/2 V 



2 + 

Then ||X"||^i is bounded by the norm of the corresponding right hand side 
avbove. By the discrete chain rule 

D\ (u-_,/^ A 0)' < 2 (u^_,f, A 0) D\u-_,/^ + At (i?>"-i/2)' • 
Estimating the first term here 

A 0) i?>;Li/2 

^) 9"-! +(""-1/2 ^0)9"+^+^" 



- I "j-1/2 A 



«^i/2V0l 



^^"-1/2 A 



^■-1/2 



AOW?" + D+P" 



This means that 



AO)!)?^!*"!!^, < C 



l£2 llg 



IP'+P" 



< c. 



Similarly, we have 



and thus 
(7.7) 



<C 
< C 



{d\u-^,,,)\c [u-^,,, V 0)' + {q-f + A 



< C. We have shown that 
{d\ [(«7_i/2A0)^ + K+i/2V0)^]} 



< C. 
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Next, using (|6.2p with f{q) ~ 0^/2 we have that 



2 



Dl [q-f = -D^ [[u-^,f, V O) [q-f] + [[u]_,/, A o) (q^) 



(7.8) 



where A" is defined in (|6.12p and 

/a.,/",, - Ax { V O) {D^q-f - (^;Vi/2 A O) [D+q^f] 

We write = F/^" + F/'" where 

y/'" = (/ - and f/^" = (/ - 6,". 

Now 1 1 F**'" 1 < 1 1 6" 1 1^1, and therefore we compute 

||AVNI,i = \A-q-\ <c\n. I Ax^ [(p;f + 



> 0. 



1/2 



1/2 



<c\m,. Ax^[|p;| + (^.;Vi/2)' 



< c. 



By (|6.13p . the norm of At{D^_^q'j)'^ is of the same order as Aa;^. Then, summing 
(j7.8p over n and j, we arrive at 



AtAx^lAxj",j <CT+ AxY^ [(. 

n,j j 

This means that 



7°)'] +o[Ax') <C. 



(7.9) 
Now let 



AtAxJ2\Yj 



< C. 



for some sequence {Kj}- G fi. Since 



we get 



Using this, 



= AxJ2\Lj\ < Ax/iC^e-"^!'-^! \Ki\ < C||i^|| 



Combining this with (|7.9p and (|7.7p we see that 
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and therefore 



/ / l^t^A.I dxdt < C. 
Jo Jr 



Hence {Pax}ax>o bounded in W^-^{{Q,T) x M). Combining this with the L°° 
estimates fomid in Lemma [4.21 yields the existence of a convergent subsequence as 
claimed in I^L^ . □ 

8. Strong convergence result 

We now show that the sequence {qax} ax>o^ P-^P - l^^s a strongly converging 
subsequence. This result is a key point of the convergence analysis. 

Lemma 8.1. Fix 1 < p < 3 and 1 < r < 1 + f • Then there exist two functions 
q e ^[^^((OjT) X R), ^2 g L[q^((0,T) X K) such that for a subsequence as Ax 0, 

(8.1) qAx^q inL^{0,T;L\R)), qAx ^ q in Ll^i{0,T) x R), 

(8.2) ql^^^ mLl,{{0,T)xR), 
for all a, 6 G M, a < b. Moreover, 

(8.3) q^{t,x) < q^{t,x) for a.e. {t,x) G (0,T) x R 
and 

(8.4) dxU = q in the sense of distributions on (0,T) x R. 
Finally, there is a positive constant C such that 

(8.5) 9(i,a:) < ^ + C||mo||^i(„), < G (0, T), x G R. 



Proof. Claims (|8.ip , (|8.2p are direct consequences of Lemmas 14.11 and 16.11 Claim 
(|8.3p is true thanks to (|8.2p and the convexity of 5 = g^, cf. Lemma [^31 while (|8.4p 
is a consequence of the definitions of g^a: and UAa:i cf. (|3.9p and (|3.10p . 

We conclude by proving (|8.5p . Fix i > 0, and let At be so small that t G /" 
with n > 0. Not that n ^ cx) as At — > 0, and t" and both tend to t. From 
the definition of gAa; and (jS.ip we have that 



qAx{t,x)^q^ + {t^t^)D\q^ 
t"+^ - t „ t-t 



97 + -r-^r' 



At ' At 

(8.6) <^f^ + cUi-^fl + C 

At v^" / V*" 

= ^+C + 2/At(t), 

where C* := C* ||uo|lifi(K) ^^.d for every t G [t", t"+^), x G Ij, with 

t - t" 1 t"+i - t 1 1 
" At t»+i At F ~ t ■ 

Observe that 

/At(i) - ^ I ^ ~ ^ ) + ^ - "t»t"+i ^ t2 ' 
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SO f'^tit) — ^-.t = V e {tn, tn+i), and in particular 



sup fAtit) = /A*(Vt"i"+l) = 

te[i",t"+i) 



< 



At 



1 



< 



Therefore (jS.Sp follows from (|8.ip and 



□ 

In view of the weak convergences stated in (IS.ip . we have that for any function 
/ G Ci(M) with /' bounded 



(8.7) 



f{qA.)^fiq) inL°°(0,T;L2(R)), 

/(<ZA.) ^ in LL((0,T)xM), l<p<3, 



where the same subsequence of Ax — > applies to any / from the specified class. 

In what follows, we let qf{q) and f'{q)q^ denote the weak limits of qAxf{qAx) 
and f'{qAx)q\,, respectively, in LL,((0,r) x R), 1 < r < |. 

Lemma 8.2. For any convex function f £ C^(R) with f bounded we have that 



(0,T)xI 



(^f{q)dtsp + uf{q)d^ip 



dx dt 



f{qoix))(p{0,x)dx 



> 



/(0,T)x 

for any nonnegative ip G C^{[0,T) x 
Proof Set 

1 



(^lr{q)q^ - qf{q) + {p- u^) /'(?)) ^dxdt, 



ipj{t) 



Ax J I 



ip{t,x)dx, (/s" 



1 



AxAt 



(p{t, x) dx dt. 



We multiply (|6.2p by AxAtip'j, sum over and take into account the convexity 
of /. After partial summations, the final result reads 

(8.9) Eo + Ei+E2+E3 + Ei + E5> 0, 

where 

Eo^AxY,fiq>l 

3 

E, ^AtAxJ2fiq7)Dr^l 

E2 =AtAxY, [(""+1/2 V O) fiq';)D+^^ + {u]_y, A o) fiq^)D_^^: 



E3=- At Ax 



"Pi 



E,^-AtAxY,A^f'{q^^)^], 

E, ^At^AxY^f" {DXq-)\l 
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By dSSD, 

||(7o - gAK,o||i2(R) ^ as Ax ^ 0, 
where qo := d^uo and qax.o dxUAx\t=o, so that 

Ax J2 /('?°)<^° ^ / /(9o(a;))v?(0, x) dx as Ax ^ 0. 

We spHt El in three parts: 
(8.10) El = Ei^i + Eia + Ei^3, 

where 



Ei^i = I f{qAx)dt(pdxdt, 

' (0,T)xR 



Due to O 



(8.11) Ei^i^ // f{q)dtipdxdt. 

J J{0,T)xS. 

Due to the boundedness of /', 

(8.12) \f{q^) - fiqAx)\ < ci(< - niDlq^l 

for each {t, x) G [i", i"^^) x Ij, where ci > is a finite constant. In view of (|6.13p . 

ii?i,2i <ciAx^lA-</'"|i^W"i / {t-ndt 

(8.13) "'J' *" 

< ciAxAi^^ |Z?r^"| < ciAt 11^^911,2 ^ 0, 

as Ax — > 0. Finally, since 

\D^(p — dtf \ < C2AX, 
for some constant C2 > 0, we have that 



(8.14) \Ei^3\<C2Ax \fiqAx)\ dxdt ^0 as Ax ^ 0. 

<j <J supp(:p) 

Clearly ([SAO]) . ([8ll|) . (l8J3l) . and (15111) imply 



(8.15) / / f{q)dtipdxdt as Ax ^ 0. 

^ J(0,T)xR 

Next, we split i?2 into four parts: 

(8.16) E2 — E2S + -E'2,2 + -^2,3 + -E'2.4, 

where 



£'2,1=// UAxf{qAx)dxfdxdt, 
'(0,T)xR 
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^2,2 - Yj 



E2,3 — ^ 



"J 



V 0) - (uAx V 0)) f{q^)D+v^ 
+ ((""-1/2 A 0) - (uAx A 0)) fiq'DD^V^ 
{UA. V 0) {f{q^) ~ fiQA.)) D+ip] 

+ (uAx A 0) (/(g^) - /((ZA.)) D^ip] dxdt, 
{uAx V 0) fiqAx) {D+(p] - dxip) 

+ (max a 0)/(qAx) {D-ifi" - dxf) 



dx dt. 



(0,T)xR 



qf[q)dx^dxdt. 



dx dt. 

Due to ((571) . 

(8.17) i;24 - 

Using the definition of wax, 

m"_1/2 - WAx < I- 
'«"_l/2 ~ WAx < ^ 

so, by Lemmas 14.11 and 17.11 

\E2,2\ < Ax^AtY (k"| + kr^l) |/(9")| + 







k"|- 









1/" i;" 






"j + 1/2 "j-1/2 



(8.18) 



+ AxAtYUl 



+ 1/2 "j-1/2 



1/(9") 11^+^" I 



< 



2Ax\\D+^\\,^ Il9ll,2||/(«)||,. 



+ 11/(9)11,. 

as Ax ^ 0. Using (|812l) and ((67T3)) . we deduce 
|£^2,3|<ciAa;^(|i^_^7| + |i?+^7 



(8.19) 



1/2 



^J-1/2 



) f {t-e)dt 



= ciAxAt'Y{\D-V'^\ + \D+v'^\) ( 
<4ciAt||i^+(^||p Hull,. ^0, 



-^i+i/2 



^j-l/2 
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as Ax 0. Finally, since \D±Lp — dt(p\ < csAx, for some constant C3 > 0, we 
obtain 



(8.20) 



\E2a\ < czAx 1 1 \qi^x\\I{q^^x)\ dtdx ^ Q asAx^O. 



/supp(ip) 

Note that (|8l6l) . ([STT]) . ([SlS]) . ((09l) . and ([QO]) imply 



(8.21) 



^2 



'(0,T)xR 

We split i?3 into three parts: 



qf{q)dx^dxdt as Ax 0. 



(8.22) 



E3 — E31 + £3,2 + E- 



3,3, 



where 



E- 



3,2 



E: 



3,3 



(0,T)xR 



(O'Ai 



E 



2 

(g"^2 



-fiQAx) - qAxfiqAx) 



ip dx dt, 



/ [qAx) + qAxfiqAx) 



(/s" c?t, 



E 



(qAx)^ „, 



/'(o'Ai;) - qAxfiqAx) 



iipj - </?) dx dt, 



Due to 



(8.23) 



£'3,1 ^ — 



II 

JJ(( 



q'f'iq) 



(0,T)xR 



— qfiq) Lp dx dt as Ax — > 0. 



Using the boundedness of /', we can estimate as follows: 



< 



2 f'iqj) - q7fi<l7) - ^^f'iqAx) + qAxfiqAx) 
|(g")^ -gLI 1^/^ „x| , (gA^iw. w 



|/'(9,")-/'(gA.)| 

\2 



^ |g,"-gA.|(|g,"| + kA.|) I (gA.)' II „|| | „ 
< 11/ IIl- ^ 11/ Hi- \1] - gAa 



AN EXPLICIT SCHEME FOR THE CAMASSA-HOLM EQUATION 



35 



Hence, taking into account (|6.13p and ()6.3 
(8.24) 



IS3.2I < A.T^ 



-TT^ 11/ = 



iirii^^ it-n\D',qj\ 



< AxAf ^ 



< At 



n,3 
1* 



Il/'ll 



\\f"\\L^ I^V?; 



wXnp (lk"ll£2 + \\qAx\\L^(t,)) 

7. 11/ IIl= 



< At 



O 



•Ax" 
At 



'Ax' 1 
Ai Aa; 



= Ax' 



as Ax ^ 0. 



Since \ip:^ - (p\ = 0{Ax) 
(8.25) 

|£;3,3| <0(Aa;) / / 



f (qAx) - QAxfiqAx) 



'supp {<p) 

We have that ([g?^ . (jOg)) . (jSTM)) . and (jS:^ imply 



dxdt ^ as Ax ^ 0. 



qf{q)dx^dxdt as Ax ^ 0. 



(8.26) E:i 

'(0,T)xR 

We spht the term into three parts: 

Ei = Ei^i + i?4_2 + i?4,3, 



where 



E4,2 



AAxf'{qAx)ipdxdt, 

'(0,T)xR 



-2^4,3 = - X] // ^AxfiqAx) {Vj - f) dxdt, 



where — Pax ~ (uax)'^- Lemmas 1 7 . 1 1 and [7^ cf. also (|8.7p . imply that 
i;4,i ^ - {P~^^) f'iq)f '^^ as Ax ^ 0. 

JJ{0,T)xS. 

Continuing, it is not hard to see that 1-^4,2 1 = 0(Ax) ^ as Ax 0. Moreover, 
since ip^ — tp — ©(Ax), E^^s ^ as Ax — > 0. Summarizing, 



(8.27) 



E4 



(P - u^) f'iq)ipdx dt as Ax 0. 
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Finally, regarding E5, due to (I6.13p . we conclude that as Ax 
(8.28) \E,\ < M^^^^,^^^^^^AfAxY,f"[q;+'^') {DXq-)\-^0. 

The lemma now follows from (US]), ([81^ . ([OT]) . ([Oe]) . ([OTI . and ([OSll . □ 

We know that {{qax)^} ^^^^ is bounded in L°°{0,T; L^{M.j) n L[oc((0'^) x I*)' 

2 

for any 1 < r < 1 + 0/2. Additionally, using (|6.2p with f{q) — we can show that 
the mapping t i— > /jjj((7A2;)^<^c?a; is equi-continuous on [0,r], for every tp e C^(R). 
Hence, in view of Lemma 12.41 



(8.29) / (qAx)'^ f dx / q^Lpdx uniformly on [0, T], 
and 

(8.30) i q^Lpdx is continuous on [0, T]. 

The statements (|8.29l) and (I8.30p hold with q^^, (7^ replaced respectively by f{qAx), 
f{qAx), for any convex fimction / G C^(]R) with /' bounded. 

Lemma 8.3. Let q and q^ be the weak limits identified in Lemma |8. 11 Then 
{qdtip + uqdx^p) dx dt + / qo{x)(p{0, x) dx 



11(0 



,.31) " 



(0,T)xl[ 



+ (p_u2^^ V(^e C:?°((0,T) x 



Proof. Starting from (j6.2p with /(g) = q, we argue as in the proof of Lemma 
to conclude the validity of (|8.3ip . □ 

The next lemma tells us that the weak limits in Lemma 18.11 satisfy the initial 
data in an appropriate sense. 

Lemma 8.4. Let q and q^ be the weak limits identified in Lemma l8. 11 Then 

\2 



lim / q {t,x)dx— / (dxUo) dx, 

(8.32) *^°Jm Jm 

lim / q'^{t,x)dx— / [dxUof' dx. 

JR JR 

Proof. The proof is similar to that in 12J. □ 

We can now wrap up the proof of the strong convergence of {'ZAa;}Aa;>o- 
Lemma 8.5. Let q and q^ be the weak limits identified in Lemma l8. 11 Then 

(8.33) Ipit, x) = q^it, x) for a.e. {t, x) £ (0, T) x M. 
Consequently, as Ax — > (along a subsequence if necessary) 

(8.34) qAx ~* q in Lf„^((0, T) x R) and a.e. in (0, T) x M. 
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Proof. By Lemma 18^ 

(8.35) + d. (um) < - i/^W + {u' - P) T{q), 

in the sense of distributions on (0, T) x R, for any convex function / e C"'^(R) with 
/' bounded Moreover, by Lemma [8?3l 



(8.36) dtq + d4uq) ^ ^q^ + - P, 

in the sense of distributions on (0,T) x M. Equipped with (lOSj) . (jOG]) . ([OS]) , 
and (18. 5|) . we can argue exactly as in Xin and Zhang [39 to arrive at (|8.33|) . In 
view of Lemma 12.31 claim ()8.34p follows immediately from ()8.33p and (17. 5|) . □ 



We now prove that the limit u satisfies (D.3) 



Lemma 8.6. For any (p £ C~((0, T) x E), 



.37) 



uifit + + P j fx dx dt = 0, 
P {(fi — (fixx) dx dt — j j [ H — (dxu)^ \ (fi dx dt 



2 

Proof. It is not difficult to establish the equation for P, since we have already 
established that dxUAx — > dxU in Li^^{{0,T) x R), cf. (|8.34p . Indeed, we have 

rT 

Pax - Vxx) dxdt 



Pax {f - fxx) dx dt 
= X! / / 'P-fxx dx dt 

+J2 [[ (-^r - ^A.) - 'Pxx) dx dt 



3-1/2 



= AtAxY,PHv]'D^D+f^l) 

+ PP [[ (((^ - + {fxx - D-D+f])) dx dt . 



Since \PAxit,x) - P^\ < CAx for {t,x) e /"„i/2, Wi ^ as Ace ^ 0. Similarly, 
since (p — ipxx is close to t^j" — D^D+cp^ in I^_i/2 Pj" is bounded, we conclude 
that W2 ^ as Ax 0. Set 

fF={<-i/2^0)\{u-^,/,V0y+^-{q-)\ 
fAx = (uAxf + ^ (dxUAxf ■ 
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Using the scheme for PJ', cf. (13.3 



/ax^ dx dt 



E 



3-1/2 



— dx dt . 



By the definition of gAa; we have that (/ak = + ^ t'^'')D\q'j for t e Hence 
SO that, assuming supp((^) C [xj^ , Xj^] for some intergers ja and jb, 



E 



(9Ax)' 1^1 dxdt 



dt 



3=]a 



= o( VAtA^ + AtAx^ 



By the Holder continuity of uax', recall that uax G C'^'^ with £ — 1 — 2/(2 + a), we 
find 



0(Ax' + At' 



and therefore -103 as Aa; 0. Hence, using (fr2|) . ([Q]) . and ([Qi]) . 



P [ip ~ ipxx) dx dt = 



lim 



-Pax (<y5 - (fxx) dxdt 



(uAx) +-^idxUAx) ]ipdxdt 



1 



{dxuf' I Lpdxdt. 



This means that the second equation in ()8.37p holds. 

To establish the first equality in (j8.37p . we derive a divergence- form version of 
the scheme p.ip . To this end, introduce the functions /v(u) = ^(u V 0)^ and 
/a(u) = ^(u A 0)^. Observe that /v and /a are piecewise C^, and the absolute 
value of the second derivatives are bounded by 1. By the discrete chain rule, 

(zi;Vi/2 V 0) D^u]^,/, = D^Mu]^,/,) + o (^Ax (i?_«;Vi/2)') 
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and 

(""+1/2 A O) D+u]^,/, = D+Mu]+,^,) + O (^Ax (i?+w;Vi/2)' 
Consequently, we can replace p.ip by 



(8.38) 



Observe that 



(8.39) 



D = 



AxD-D+=D+~D., /v + /a = 



Using these identities, we can restate ()8.38p as 



(8.40) 



^+""+1/2 + D- 



'j + l/2 



^■ + 1/2 



4 

O ( Ax 



+ ^(/vKVi/2)-/a(";Vi/2)) 

+ ^(./a(";Vi/2)-/v(^;Vv2) 



;"+i/2 



2 / \ 2 

D+ul 



''3 + 1/2 



Using, cf. (jOg)) . 

(/vl^i^'+l/^) - /a«+1/2)) + ^+ (./a (""+1/2) - /v(^i "+1/2)) 

= AxD^D+f^{u';^„^) - AxD^D+Mu^^,/^), 
equation l|8.40p becomes 

^(-"+1/2)' 



+ S-+1/2 + ^ 



(8.41) 



V 



A: 



Now fix G C^((0,T) X R) and define ip'^ as before, cf. (HJ]). Multiplying ((OTI) 
by (p^ AtAx and performing partial summations gives 



AtAx^u'^^y^D\ip'^+AtAxY^ ^ ^'"l''^' Dip]+AtAxJ2P"D^ip] = 0{Ax), 



El 



E3 
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by using (|4.2p . We have UAa 
this and (|7.2p . we compute as follows: 



^i-l/2 



< C (^Ax\q:^\+ At 



n* 1/" 

+ i+i/2 



dt 



Using 



^1 ^ J J '^Ax^t dxdt + '^ J J (w"+i/2 ~ ■"Ax) t/5t dx 

+ E ""+1/2 / /„ (^+^" - ^0 dt 

71, j 

= / UAxftdxdt + 0{Ax) / / uiptdxdt as Ax ^ 0. 



In the same way, equipped with (|7.2p and (|7.6|1 . we can show that 



-ifxdxdt, i?3 



P(px dx dt as Ax — > 0. 



thus proving the first equality in (|8.37p . This concludes the proof of the lemma. □ 

9. Numerical examples 

We have tried the difference method presented here on several examples, and 
in doing this found that the first order method analyzed in this paper exhibits 
very slow convergence, and thus requires a very small mesh size Aa; to compute 
reasonable solutions. This is not surprising and appears to be the case with other 
schemes in the literature as well. Therefore we have implemented a second order 
extension of the method. This second order extension is based on the conservative 
version of the scheme 

(9.1) i?>"+l/2 + D- [{u]+,/^ V O) ^7+1/2 + A O) u]+y^ + f^+i" 

= D- (uJ+1/2 V O) <_i/2 + (,.;Vl/2 A O) 

which can be viewed as a balance equation with a flux across x = Xj_^_i/2 given by 



^.7 + 1 



i+1/2 



'■j+1/2 



vou 



'■j+1/2 



AOU 



^J+3/2 



Taking this viewpoint, we define the second order finite volume scheme by 



D^F 



n+l/2 



(9.2) , ^-.,+1/2 - 

Here u^^^j2 ^ &rst order approximation of the value at the point x — Xjj^_ii2, 
i = + At/2. This approximation is found by setting 

and then using the scheme p.ip for half a time step (i.e.. At/ 2). This scheme 
is a formally second order accurate finite volume approximation, and this simple 
adaptation produces significantly more accurate approximations. 

In Figure [1] we show the approximations calculated by the first order scheme 
()3.ip and the second order scheme (|9.2p for the single peakon example. In this case 
the exact solution reads 



u(x, t) 
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Figure [T] shows the solutions calculated using 2^ equally spaced grid points in the 
interval [—10,30] for t ~ 20. We see that the second order method is much more 
accurate than the first order method. In passing, we note that we have not used 



Single peakon approximation at f=20 




-5 5 10 15 20 



Figure 1. Approximations using Aa: = 40/2^, at i = 20 to the 
single peakon. 

the strict CFL-condition (|4.ip . but the more natural condition 

At < max Ax. 

3 

This holds for the second order scheme as well. 

In order to investigate the convergence properties of the two methods, we com- 
puted errors in for the two schemes. Table [T] shows the computed errors in 
the case of a single peakon under mesh refinement. In this context, the error is 
defined as 

error ^ Ax'^ \uAxixi,t) - u{xi,t) \ , 

i 

where u is the exact solution. We used t = 20 and Aa; — 40/2*^ for A; = 5, ... , 12. 
As expected, and as reported in [12j , the first order method converges very slowly. 
One other notable feature of Table [T] is that the second order method seems to 



k 


5 


6 


7 


8 


9 


10 


11 


12 


13 


2nd 


2.92 
5.36 


3.23 
5.17 


3.41 
3.29 


3.53 
1.27 


3.57 
0.60 


3.51 
0.36 


3.32 
0.21 


3.01 
0.13 


2.64 
0.09 



Table 1. errors for the single peakon case, dX t — 20, for 
X e [-10, 30], Aa; = 40/2*^, fc = 5, . . . , 13 



converge at a rate slightly less than 1. 

The two-peakon solution is considerably more complicated than the single peakon, 
and this is also a much harder challenge computationally, see e.g., [T] and [33]. We 
use the two-peakon solution given by 

(9.3) u(a;,t) = mi(t)e-l^-^i(*)l + 7n2(t)e~l="-^^(*)l , 
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with 

Mt) = log {J^^) , Mt) - log (40e-- + 60e(--)/^) 
e(t-io)/2 + 6 e(t-io)/2+2 
"^ift)- 2e(«-i°)/^ + 3 ' e(*-i")/^ + 3 - 

These formulas were taken from [35j . Figure [2] shows a contour plot of the approxi- 
mate solutions found by using the first and second order methods, and Ax = 40/2^*^ 
for X £ [—15,25] and t E [0,25]. We see that the interaction between the two 




Figure 2. Approximations to (|9.3p using Ax = 40/2^". Left: first 
order method( |3.ip . Right : second order method (|9.2|) . 

peakons is poorly represented by the first order method. Both the location as well 
as the magnitude of the peaks are far from the correct value. This is also illus- 
trated Figure [3] where we show the approximations using Aa; = 40/2^ at t — 25. 
We have also calculated errors for the two-peakon case. Indeed for Aa; > 40/2^^, 



Two peakon approximalion at f=25 




Figure 3. The approximations to at i = 25 and Aa; = 40/2^. 
the first order method did not seem to converge, and in order to give meaningful 



AN EXPLICIT SCHEME FOR THE CAMASSA-HOLM EQUATION 



43 



k 


8 


9 


10 


11 


12 


13 


2nd 


4.56 
1.88 


3.64 
1.04 


3.97 
0.63 


4.18 
0.38 


4.05 
0.22 


3.70 
0.16 



Table 2. errors for the approximation to l|9.2p . t = 25, x G 
[-15, 25], Ax = 40/2^ A: = 8, . . . , 13. 



answers, this method demands very fine discretizations. These results are reported 
in Table [21 For Ax > 40/2* none of the methods gave satisfactory results. 

In our final example we choose initial data corresponding to a peakon-antipeakon 
collision: 

(9.4) uoix) = - tanh(6) (e^l^+^^^^l - e^l^-''^^)!) , 

where y(t) — log(cosh(i)). In this case we have a "peakon anti-peakon collision" 
at t = 6. In Figure 2] we exhibit the approximations generated by the first order 
(left) and the second order method for t G [0,10] and Aa; = 24/2^^. It is clear 
that the first order scheme generates the dissipative solution, and for t larger than 
the collision time, the first order approximation vanishes. Regarding the second 
order approximation, it seems to continue as a peakon moving to the right, and an 
anti-peakon moving to the left. The magnitudes and speeds of these features are 
however far from the conservative solution, and we have indicated the conservative 
solution in the right hand figure. 




Figure 4. The numerical solutions to the initial value problem 
(|9.4p . Left: first order method, right: second order version. 
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